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SECTION  I 
INTRODUCTION 

When  a  vehicle  reenters  the  earth  atmosphere,  it 
passes  through  the  different  flow  regions  from  free  mole¬ 
cule  to  continuum.  This  has  created  many  new  and  difficult 
problems  for  engineers  in  many  fields.  As  far  as  heat 
transfer  is  concerned,  the  most  important  problem  is  that 
of  excessive  temperature,  due  to  skin  friction,  attained 
at  the  surface  of  the  solid  body  traveling  at  extremely 
high  velocity. 

Heat  transfer  has  been  of  interest  to  the  engineers 
and  rocket  sled  designers  for  many  years.  With  the  advent 
of  reentry  vehicles,  velocities  of  rocket  sleds  have  had 
to  be  increased  to  maintain  their  usefulness  as  a  suitable 
testbed.  Therefore,  heat  transfer  and  temperature  dis¬ 
tribution  have  become  the  primary  problem  rather  than  just 
of  interest  to  sled  designers.  Three  independent  primary 
heat  sources  affect  hypersonic  sleds:  heat  generated  by 
slipper-rail  friction;  heat  generated  by  oxidation  of  sled 
components;  and  aerodynamic  heating.  This  last  effect 
will  be  analyzed  in  this  thesis.  The  temperature  distribu¬ 
tion  is  necessary  in  the  design  of  rocket  sleds  from  a 
structural  standpoint  due  to  thermal  stresses,  as  well  as 
temperature-time  profiles  in  the  sled  for  temperature- 
sensitive  instruments  and  components. 
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Since  the  sled  trajectory  is  in  the  atmosphere  near 
the  earth,  the  mean  free  path  of  gas  molecules  is  rela¬ 
tively  very  short  in  comparison  with  the  actual  dimensions 
of  a  solid  body  in  movement.  In  conventional  aerodynamics 
and  heat  transfer,  air  is  usually  treated  as  a  homogeneous 
continuous  medium.  That  is  to  say  the  gas  has  a  definite 
density,  pressure,  and  temperature. 

The  object  of  this  program  is  to  establish  a  method 
of  predicting  boundary-layer  conditions,  the  heatflux, 
and  the  time-temperature  distribution  in  the  rocket  sled 
skin.  The  thermal  model  examined  in  this  thesis  assumes 
thick-skin  solution  where  the  skin  is  of  any  composite 
structure  made  of  discrete  layers  of  material  whose  prop¬ 
erties  are  temperature  dependent  and  vary  from  layer  to 
layer . 

The  boundary-layer  conditions  are  obtained  from  the 

flat-plate  solution  modified  for  conical  flow  by  the 

Mangier  transformation^  and  for  compressibility  by 

(2 ) 

Eckert's  reference  temperature  method.  These  boundary- 
layer  properties  are  used  to  determine  the  heatflux  and 
temperature  profiles  in  wedges  and  protuberances  down¬ 
stream  on  the  sled. 

This  program  is  specific  in  that  the  equations  have 
been  applied  to  a  blunt-nosed  conical  vehicle  with  the 
starting  point  at  ignition  using  ambient  conditions.  With 
some  modifications,  the  user  could  easily  adapt  the  pro¬ 
gram  to  analyze  any  sled  configuration. 


SECTION  II 


THEORY 


1.  General 

This  computer  program  is  a  self-contained  program  in 
that  it  can  generate  the  required  aerodynamic,  heat  trans¬ 
fer,  and  trajectory  data  without  any  additional  prelimi¬ 
nary  calculations.  The  program  calculates  loads,  heat 
flux,  temperature-profile,  and  boundary-layer  properties 
as  well  as  thrust,  velocity,  and  displacement  all  as  a 
function  of  time. 

At  the  initiation  of  a  sled  run,  the  vehicle  is  lo¬ 
cated  on  the  track  at  any  station  with  the  desired  ambient 
atmospheric  conditions  (i.e.,  temperature,  pressure). 

That  portion  of  the  sled  under  study  is  divided  into  a 
number  of  segments  or  stations,  in  this  case  the  cone  or 
wedge,  and  cylinder.  The  heating  rate  at  each  station  is 
computed,  and  from  this  the  temperature  profile  is  obtained. 
In  addition  to  atmospheric  conditions  required,  the  fore 
and  aft  radii,  slant  length,  and  transition  Reynolds  num¬ 
ber  of  each  body  segment  must  be  supplied.  For  each  inte¬ 
gration  step,  At,  a  complete  new  set  of  data  is  computed 
for  each  segment. 

2.  Mathematical  Model 

Before  proceeding  with  the  theoretical  development, 
it  may  be  worthwhile  to  digress  briefly  and  discuss  types 
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c.  HYPERSONIC  MONORAIL  IMPACT  SLED  (20  inch  DIA) 


d  HYPERSONIC  MONORAIL  SLEO  FOR  TESTING  AIR  BREATHING  PROPULSION  SPECIMENS 


Figure  1.  Examples  of  Hypersonic  Monorail 
Sled  Configurations 


of  Sieds  used  in  this  study  and  the  mathematical  model. 
Referring  to  Figure  1,  there  are  five  typical  examples  of 
hypersonic  monorail  configurations.  Examination  of  these 
sleds  illustrates  those  components  of  the  sled  or  payloads 
most  vulnerable  to  aerodynamic  heating  which  are  subjected 
to  stagnation  conditions,  and  those  areas  subjected  to 
multiple  shock  interaction.  Special  trouble  spots  are 
body  transitions  from  cones  to  wedges,  or  cylindrical  body 
shapes . 

This  study  is  concerned  with  the  stagnation  temper¬ 
atures  of  a  cone  as  shown  in  Figure  2.  Also  of  concern 
are  temperature  profiles  at  various  stations  along  the 
cone . 


Figure  2.  Sled  Configuration  Used  in  this  Study 

The  nose  cone  (Figure  3)  is  divided  into  N  number  of 
segments  or  truncated  cones  with  the  temperature  calculated 
at  the  midpoint  and  assumed  constant  for  that  segment. 
One-dimensional  heat  flow  is  used  since  the  time  of  flight 
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STA  STA  STA  STA 


is  relatively  short  and  the  heat  flow  from  the  boundary 
layer  is  much  greater  than  heat  flow  along  the  skin. 

3.  Aerodynamic  Heating 

Hypersonic  vehicles  in  general  may  be  divided  into 
two  categories:  The  first  category  comprises  those  ve¬ 
hicles  which  must  be  propelled  through  the  atmosphere  to 
sustain  flight.  The  second  category  encompasses  those 
that  have  a  vast  store  of  kinetic  and  potential  energy 
which  must  be  dissipated  during  their  flight  through  the 
atmosphere.  In  the  former  category  is  the  aerodynamic 
vehicle  and  in  which  the  rocket  sled  belongs.  In  either 
case,  some  of  this  energy  is  expended  against  drag  forces 
which  assume  the  form  of  skin  friction  or  pressure  drag, 
and  produce  heat. 
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Heat  Transfer  Coefficient 

To  determine  the  heat  transfer  coefficient,  h  ,  the 

c 

nature  of  the  flow  can  be  described  by  the  Reynolds  number, 
which  is  a  dimensionless  measure  of  the  ratio  of  inertial 
to  viscous  forces.  For  this  study,  the  flow  was  assumed 
turbulent  from  first  movement,  because  of  the  extremely 
rapid  acceleration  of  the  sled  and  from  the  ground  inter¬ 
ference  effects.  The  ground  interference  phenomena  asso¬ 
ciated  with  rocket  sleds  are  due  to  their  moving  at  high 
speeds  close  to,  and  parallel  with,  the  ground  plane. 

This  so-called  ground  interference  phenomenon  will  be  evi¬ 
denced  by  a  change  in  the  flow  velocity  and  pressure  fields, 
in  the  region  between  the  sled  and  the  ground,  when  com¬ 
pared  to  the  flow  and  pressure  fields  existing  in  the  ab¬ 
sence  of  a  ground  plane. 

In  high  speed  flow,  at  least  two  additional  param¬ 
eters  must  be  considered.  These  are  Mach  number  and 
Prandtl  number.  The  parameter  Mach  number  describes  the 
influence  of  compressibility  on  heat  transfer  and  flow 
phenomena.  The  Prandtl  number  is  defined  as  the  ratio  of 
heat  storage  to  heat  conduction  of  a  gas.  The  Reynolds, 
Prandtl,  and  Mach  numbers  are  the  parameters  governing 
convective-heat  transfer  as  long  as  the  flow  region  may 
be  treated  as  a  continuum. 

Therefore,  the  variable  affecting  heat  transfer  by 
turbulent  forced  convection  to  or  from  a  gas  flowing  over 
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a  flat  plate  or  cone  can  be  obtained  from  the  analogy*  be¬ 
tween  momentum  and  heat  transfer.  This  analogy  was  first 
developed  by  Osbourne  Reynolds. 

The  following  equation  is  known  as  the  Reynolds  anal¬ 
ogy  between  momentum  and  heat  transfer,  or  between  the 
fields  of  velocity  and  temperature 


a 

A 


TCpa  dT 
v  "  X 


(1) 


where  the  Prandtl  number  is  not  equal  to  unity.  (See  Ap¬ 
pendix  A. ) 


v  =  kinematic  viscosity  -  u/P 
a  =  thermal  diffusivity  -  k/pCp 

The  basic  analogy  was  later  improved  by  Prandtl ^ 
and  others,  whose  additional  refinements  were  particularly 
applicable  to  forced  convection  flow  over  flat  plates. 
Colburn ^  established  the  following  correlation  between 
heat  and  momentum  transfer  based  upon  Reynolds'  analogy. 

h,v  =  0.0296  ~  (P  ) 1/3  (R  )  °* 8  (2) 

CX  X  i  6 

which  is  equation  (A-12)  from  Appendix  A. 


•See  Appendix  A. 
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The  solution  of  this  equation  for  the  heat  transfer 
coefficient  in  conical  flow  is  obtained  by  equation  (A-13) 


cx 


»  0.0339  £ 


(Rex)°-8(Pr) 


1/3 /u* 


.  * ,0.2  ,  *  ,0.8 

e  e 


Good  correlation  with  experimental  data  is  obtained 
if  the  gas  properties*  are  evaluated  at  the  reference 
temperature  and  the  velocity  is  taken  as  the  boundary- 
layer  edge*  velocity. 

The  surface  of  the  vehicle  may  also  receive  heat  by 
radiation  from  the  hct  gas  in  the  shock  layer  in  addition 
to  that  produced  by  convection.  Since  the  sled  components 
considered  in  this  study  were  polished,  the  rate  of  radi¬ 
ation  heat  transfer  is  negJ  j.gible  compared  with  that  of 
aerodynamic  convection  and  can  be  ignored. 

Heat  Flux  Equation 

Experiments  in  high-speed  flow  have  verified  that 
the  magnitude  and  direction  of  heat  flow  at  the  surface 
do  not  depend  on  the  difference  between  the  wall  tempera¬ 
ture  and  the  free-stream  temperature  as  in  low-speed  flow, 
but  rather  on  the  wall  temperature  and  the  adiabatic  wall 
temperature.  It  is  apparent  that  the  determination  of 
the  adiabatic  wall  temperature  will  be  of  prime  importance 
in  the  calculation  of  heat  transfer,  since  the  transference 

*See  Appendix  B. 
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of  heat  to  or  from  the  wall  will  depend  upon  whether  the 
skin  temperature  is  above  or  below  the  adiabatic  wall 
temperature.  The  adiabatic  wall  temperature  can  conven¬ 
iently  be  expressed  in  terms  of  a  dynamic  temperature  rise 


where  the  gas  properties  are  evaluated  at  boundary-layer 
edge  conditions.* 

Reference  7  shows  that  for  practical  purposes  the 
recovery  factor  for  turbulent  flow  is 

RF  =  /FT  (5) 


which  is  based  upon  the  assumption  of  con' tant  eddy  dif- 
fusivity  for  momentum  and  heat. 

Therefore,  the  unit  surface  convective  heat  rate  for 
high-speed  flow^  is 


h  (T... 
cx  AW 


V 


(6) 


As  noted  in  Appendix  A,  there  is  always  a  laminar 
sublayer  and  a  buffer  zone  under  the  turbulent  layer.  The 
rate  of  momentum  exchange  varies  from  layer  to  layer.  In 
the  laminar  sublayer  near  the  vehicle  wall,  momentum 


*See  Appendix  B. 


exchange  is  negligible,  and  heat  flows  almost  solely  by 
molecular  conduction;  in  regions  sufficiently  far  from  the 
interface  and  outside  the  laminar  layer,  momentum  exchange 
is  the  prevailing  process  of  heat  transmission.  Therefore, 
in  using  Reynolds'  analogy,  the  laminar  sublayer  has  been 
neglected . 

4.  Conductive  Heat  Transfer 

The  numerical  computation  and  solution  of  the  Fourier 
conduction  equation  for  unsteady-state  conditions  is,  in 
most  cases,  rather  complex  which  involves  solving  a  trans¬ 
cendental  equation.  This  solution  is  obtained  by  deter¬ 
mining  an  infinite  number  of  eigenvalues  or  characteris¬ 
tic  values  from  a  plot  of  curves.  Due  to  this  complexity, 
it  is  more  convenient  to  solve  unsteady-state  problems 
graphically. 

/  8 ) 

A  method  was  developed  by  Binder  in  1911,  and 

IQ) 

later,  in  1924,  it  was  improved  by  Schmidt.  The  prin¬ 
ciple  of  finite  difference  is  applied  to  the  solution  of 
the  Fourier  conduction  equation  (a  partial  differential 
equation  that  is  second  order  in  space  and  first  order  in 
time) 


3T 

at 


(7) 


for  unsteady-i^tate  heat  transfer. 


r 


If  the  finite  differences  AT,  Ax,  and  At  are  used 
instead  of  3T,  3x,  and  3t,  the  general  equation  for 
unsteady-state  heat  conduction  becomes 


A.  T  A2T 

t  x 

w~  -  a  77~i 

(Ax) 


(8) 


where  A  T  and  A  T  are  the  finite  increments  of  tempera- 
ture  with  respect  to  the  time  t  and  the  distance  x. 

A  means  of  applying  the  Schmidt  graphical  method  for 
solving  an  unsteady  heat  conduction  is  presented  in  Ap¬ 
pendix  C.  A  thorough  discussion  of  the  theory  and  proof 
will  be  found.  The  following  equation 


T  -  T 

n,t+l  n,t 


aAt 

(Ax) 


*-(T  ,  i  ,  —  2T  .  +  T  .  .  )  (9) 

2  n+l,t  n,t  n-l,t 


is  the  transform  of  equation  (8)  to  the  numerical  solu¬ 
tion.  This  method  Is  widely  used  because  it  gives  an  it¬ 
erative  profile  of  the  temperature  change.  Also,  it  is 
quite  flexible  in  that  difficult  boundary  conditions  can 
be  handled  easily.  The  temperature  throughout  a  wall  or 
slab  can  now  be  computed  for  any  later  time  if  the  initial 
distribution  is  known. 


Vehicle  Therm?!  Model 

To  determine  the  time-temperature  profile  in  the  sled 
skin,  the  Fourier  conduction  equation  was  written  as  a 
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finite  difference  equation  and  solved  numerically.  As 
mentioned  in  Appendix  C,  the  wall  is  divided  into  a  num¬ 
ber  of  lamina  with  known  thickness,  specific  heat,  and 
thermal  conductivity,  and  with  a  known  initial  temperature 
distribution.  Figure  4  illustrates  the  model  with  its  as¬ 
sociated  nomenclature. 


TW  (!) 
TW  (2) 
TW  (3) 


AIRFLOW 


TW  (N) 

TW  (N  +  l) 


¥ 


Figure  4.  Model  for  Wall  Temperature  Profile 


The  difference  equation  (9)  for  computing  the  temper¬ 
ature  at  position  n  for  time  t  +  1  is 


T  =  t  + 

n,t+l  n,t 


aAt 


(Ax) 


2(Tn+l,t  -  2Tn,t  +  Vl,t>  (10) 


in  this  form,  the  material  properties  cannot  vary  from 
layer  to  layer,  the  above  equation  must  be  written  in  the 
form 
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T  =  T  .  +  At 

n,t+l  n,t 


n-1 


(Ax  , ) 
n-1 


t(T  _  -  T  ) 

2  n-1 ,  t  n  t 


n 


(Ax  ) 
n 


"5^Tn,t  Tn+l,t^ 


(11) 


where  a,  the  thermal  diffusivity  may  be  replaced  by  its 
defined  equivalent 


a  = 


(12) 


Equation  (11)  allows  certain  terms  to  be  grouped  together 
as  the  heat  is  conducted  through  each  lamina.  With  this 
form  of  grouping,  the  graphical  method  can  easily  be  ex¬ 
tended  to  permit  the  convective  heat  flux  to  be  used  in 
solving  for  the  wall  temperature.  The  boundary  conditions 
for  this  type  of  problem  can  be  stated:  at  any  instant 
of  time,  the  convective  heat  flux  flowing  from  the  fluid 
to  the  wall  must  be  equal  to  the  heat  flowing  by  conduc¬ 
tion  from  the  surface  toward  the  first  nodal  interface. 
This  can  be  expressed  symbolically  as 


x-o 


hcx(TAW  "  TWJ  ~k  3x 


x=o 


(13) 


In  the  graphical  solution,  the  boundary  condition  ex¬ 
pressed  by  equation  (13)  transformed  into  a  difference 
equation  as 
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/  AT  \fc  _  hcx(TAK  "  TW  ) 

Ux/o  "  -k 


or  expanded  into 


(14) 


Tn-l,t+l  Tn-l,t  +  (pc  Ax)n_^  [^cx  ^  TAW  Tn-l,t^] 


2Atk 


(pc  Ax2) 


TT— <Vl,t  -  T„,t» 


(15) 


n-1 


The  backface  boundary  condition  can  be  specified  in 
several  ways.  If  internal  heating  (or  cooling)  is  pres¬ 
ent,  the  backface  may  be  held  constant  or  varied  in  a 
specified  manner.  A  conservative  assumption  is  that  no 
heat  flows  through  the  last  lamina.  This  assumption  will 
cause  somewhat  higher  temperatures  than  would  be  encoun¬ 
tered  with  a  cooled  backface  or  other  heat  sink. 


5.  Stagnation  Heating 

A  knowledge  of  the  local  heat  transfer  at  the  for¬ 
ward  stagnation  point  of  blunt-nosed  bodies  moving  through 
the  atmosphere  is  of  considerable  importance  because  of 
the  high  rate  of  heat  transfer  which  may  occur.  In  this 
report,  an  analysis  is  made  for  the  unsteady,  forced  con¬ 
vection  using  the  modified  Lee's  equation  (Ref.  10)  for 
an  axisymmetrical  blunt-nose. 
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This  modification  was  developed  by  Fay  and  Riddell/^ 
They  also  numerically  solved  (by  successive  approximation) 
the  transformed  momentum  and  energy  boundary  layer  equa¬ 
tions  for  the  axisymmetric  stagnation  point.  In  perform¬ 
ing  the  numerical  iterations,  the  following  assumptions 
were  made: 


1)  Lewis  number  =  1 

2)  Pr  =  0.71 

w 

3)  y  =  1.10-1.20  at  high  temperature  which  is  the 
mean  ratio  of  specific  heats  behind  the  bow  shock  wave 

4)  pu  =  peue 

For  axisymmetric  bodies  which  are  not  too  blunt,  the 
velocity  gradient  at  the  stagnation  point  can  be  obtained 
from  the  following  modified  Newtonian  flow  equation: 


(16) 


as 


(v- 

'  oo 


dve 

3e~ 


J_\1/2 

M2  f 

CO  f 


(17) 
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where 


For  very  blunt  bodies  experimental  pressure  distributions 

must  be  used  to  determine  the  velocity  gradient.  Boison 
(12) 

and  Curtiss  have  correlated  experimental  stagnation 
point  velocity  gradient  measurements  with  equation  (17) 
for  a  range  of  shapes. 


Figure  5.  Blunt-Nose  Flow 

Substituting  equations  (16)  and  (18)  into  Lee's  equa¬ 
tion  (Ref.  10) 


q  =  0.94/p  p  vAT”  •  h  •  G  (M  ,  •>  ,  y  ) 
ws  ws  00  se  00  1  » 

s 


(19) 
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where 


h  =  V2/2g  J 
se  00  3c 

=  1.4  ratio  of  specific  heat  for  air 

p  =  density  evaluated  at  stagnation  reference  con¬ 
ditions 

U  =  viscosity  evaluated  at  stagnation  reference 
ws  conditions 


Therefore,  equation  (19)  can  be  used  as  a  reasonable 
approximation  for  the  stagnation-point  heat  transfer  for 
any  rounded-nose  blunt  body. 
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SECTION  III 


EXPERIMENT 

1.  Test  Description 

Three  sled  tests  were  instrumented  to  obtain  temper¬ 
ature  data  for  verification  of  the  analytical  method  pre¬ 
sented  in  this  report.  The  instrumentation  consisted  both 
of  active  and  passive  systems. 


Figure  6.  Probe  Installation  for  First  Sled  Test 
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Due  to  payload  space  and  telemetry  channel  availa¬ 
bility,  only  two  one-stage  rocket  sleds  were  instrumented 
with  the  thermocouple  probes.  The  first  sled  test  had 
two  temperature  probes  installed  as  shown  in  figure  6. 

The  second  sled  test  had  one  probe  on  the  cone  and  a 
transducer  for  heat  flux  measurement.  See  figure  7  for 
instrumentation  location. 


Figure  7.  Probe  and  Gauge  Installation  for 
Second  Sled  Test 

The  third  sled  was  run  on  a  small  two-stage,  recov¬ 
erable  rocket  sled.  Since  telemetry  was  not  available,  a 
passive  temperature  device  or  Temp-Plate*  was  used. 


^William  Wahl  Corp.,  Santa  Monica,  Calif. 
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2.  Instrumentation 

The  temperature  measuring  device  used  with  the  active 
system  was  the  MO-RE “  surface  temperature  probe.*  This 
probe  has  been  used  for  rapidly  varying  surface  tempera¬ 
tures  such  as  encountered  in  combustion  chambers,  gun 
barrels,  etc.  The  device  consists  of  a  coaxial  center 
lead  wire  housed  in  a  thick-walled  metal  tube  separated 
by  an  0.0005-inch  aluminum  oxide  insulation.  The  aluminum 
oxide  is  effective  to  2000°F  for  steady-state  tempera¬ 
tures,  and  over  3000°F  for  transient  temperatures.  The 
thermal  junction  is  formed  by  a  vacuum-deposited  metallic 
plating.  Two  and  one-half  foot  lead  wires  are  welded  to 
the  coaxial  probe  and  embedded  in  a  refractory  insulating 
cement.  To  provide  shielding  and  resistance  to  rough 
handling,  the  lead  wires  are  covered  with  wire  overbraid. 

The  thermocouples  are  provided  with  the  same  met3l 
as  the  probe  tube  into  which  they  are  placed.  The  sta¬ 
tion  location  and  installation  are  shown  in  figure  6  for 
the  first  sled  test. 

The  heat  flux  transducer  is  a  Gardon  gauge’  which 
provides  a  self-generated  millivolt  output  in  direct  pro¬ 
portion  to  the  thermal  energy  absorbed  by  the  sensor. 

The  sensor  is  a  olackened  thin  disk  of  low  thermoconduc¬ 
tivity  material.  The  disk  is  connected  at  its  edges  to 

*Heat  Technology  Laboratory,  Inc.,  Huntsville,  Ala.  35805. 
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a  larger  mass  of  material  (heat  sink)  having  a  thermoelec 
trie  potential  markedly  different  from  that  of  the  disk 
material.  One  end  of  a  fine  wire  of  heat  sink  material 
is  butt-welded  to  the  center  of  the  disk;  another  wire  of 
heat  sink  material  is  connected  to  the  larger  mass.  This 
method  forms  a  differential  thermocouple  between  the  cen¬ 
ter  and  edge  of  the  disk  which  is  illustrated  in  figure  8 

LOW  CONDUCTIVITY 


Figure  8.  Keat  Flux  Transducer 

When  the  disk  is  exposed  to  heat  flux,  the  heat  ab¬ 
sorbed  by  the  disk  is  transferred  radially  to  the  heat 
sink.  An  equilibrium  temperature  which  is  proportional 
to  the  energy  absorbed  is  rapidly  established  between  the 
center  and  edge  of  the  disk.  The  heat  flux,  q,  is  evalu¬ 
ated  by 

q  =  K  •  E  (20) 
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where  K  is  a  constant  determined  from  calibration  and  E 
is  the  equilibrium  electric  potential  between  the  two 
lead  wires  in  the  instrument  output. 

The  passive  system  for  obtaining  the  temperature  was 
a  small  temperature-indicating  sticker  called  a  Temp- 
Plate.  The  Temp-Plate  contains  heat-sensitive  elements 
that  are  hermetically  sealed  in  laminated,  high-tempera¬ 
ture-resistant  plastic.  These  stickers  are  bonded  to  any 
surface  for  which  the  temperature  record  is  desired. 
Temp-Plates  will  indicate  calibrated  temperatures  within 
an  accuracy  of  +1  percent.  Stated  temperatures  are  indi¬ 
cated  by  a  change  of  the  indicator  from  pastel  color  to 
solid  black  (see  figure  9)  .  The  change  to  black  is  irre¬ 
versible  and  cannot  be  altered,  serving  as  a  positive 
record  of  the  temperature  exposure.  Table  I  illustrates 
temperature  selection  range  used  for  Temp-Plates. 


TABLE  I 


TEMPERATURE  SELECTION  RANGE 

Temp-Plate  #1  500-600-700-900-1000-1100°F 

Temp-Plate  #2  500-600-700-750-  800-1100°F 

These  selection  ranges  allowed  the  expected  temperature 
to  be  covered  within  ±10°F.  Although  these  Temp-Plates 
will  not  give  temperature  profiles  as  a  function  of  time, 
they  will  give  maximum  temperatures. 


SECTION  IV 


COMPARISON  OF  EXPERIMENTAL 
AND  THEORETICAL  DATA 


1.  Results 

This  section  presents  the  comparison  of  the  experi¬ 
mental  data,  obtained  from  three  sled  runs,  with  the  theo¬ 
retical  predictions  (see  Table  II)  .  The  results  of  this 
study  are  presented  graphically  in  figures  10  through  18. 
Predictions  from  the  computer  program  are  reasonably 
close;  also  it  shows  consistently  the  same  general  trend 
as  the  experimental  data.  The  only  exception  to  this  is 
that  the  initial  rise  of  the  theoretical  prediction  does 
not  agree  with  the  experimental.  This  is  because  the 
boundary- layer  edge  conditions  are  considered  to  be  ambi¬ 
ent  below  Mach  1.  The  equations  used  to  develop  sharp- 
cone  boundary-layer  edge  conditions  are  only  valid  for 
Mach  numbers  greater  than  one. 

To  substantiate  the  validity  of  the  sharp-cone 
boundary- layer  edge  conditions  used  in  this  study,  a  com¬ 
parison  of  this  method  developed  in  reference  3  and  flight 
data  from  tne  Fate  reentry  vehicle  tests  are  presented  in 
reference  14.  The  results  of  this  comparison  show  excep¬ 
tionally  good  agreement  on  the  straight  cone.  When  the 
results  of  this  report  and  reference  14  are  observed,  de¬ 
viations  are  apparent  between  theoretical  results  and  the 
measured  data  at  the  beginning  of  the  cone,  behind  the 
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TABLE  II 
TEST  CONDITIONS 

Test  Run  March  1968 

Ambient  pressure  -  25.955  in.  Hg 
Ambient  temperature  -  44 °F 

Maximum  velocity  -  Actual  V  =  4516  fps 

Predicted  V  =  4487  fps 

Test  Run  May  1969 

Ambient  pressure  -  25.755  in.  Hg 
Ambient  temperature  -  74°F 

Maximum  velocity  -  Actual  V  =  4507  fps 

Predicted  V  =  4494  fps 

Test  Run  September  1968 

Ambient  pressure  -  25.700  in.  Hg 
Ambient  temperature  -  81°F 

Maximum  velocity  -  Actual  V  =  5840  fps 

Predicted  V  =  5772  fps 

stagnation  region.  This  is  most  likely  due  to  unestab¬ 
lished  turbulent  flow  and  bluntness  effects,  which  re¬ 
sulted  in  low  heat  transfer.  The  da  ;a  on  the  remainder 
of  the  models,  however,  provide  a  good  comparison  with 
theoretical  predictions. 

Flight  data  were  not  available  for  the  stagnation 
point;  therefore,  wind-tunnel  data*  were  used  exclusively 

‘Unpublished  data  supplied  by  Arnold  Engineering  Develop¬ 
ment  Center,  Von  Karman  Facility,  Arnold  A.F.S.,  Tenn. 
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Figure  10.  Temperature  Profile,  Sta.  1  Ft.,  March  1968 


Figure  11.  Temperature  Profile,  Sta.  2  Ft.,  March  1968 
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Figure  12.  Temperature  Profile,  Sta.  1  Ft.,  May  1969 


THEORETICAL 

EXPERIMENTAL 


Figure  13.  Heat  Flux  vs.  Time,  Sta.  2  Ft.,  May  196:# 
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Figure  15.  Velocity  vs.  Time,  Fate  Flight 
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Figure  16.  Temperature  Profile,  Fate  Flight 
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Figure  17.  Temperature  Profile,  Fate  Flight 
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Figure  18.  Stagnation  Point  Temperature 
for  Hemisphere-Cylinder 
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for  this  comparison.  The  comparison  using  the  wind-tunnel 
data  presented  in  figure  18  is  for  a  sphere-cone  configu¬ 
ration  in  air  at  a  free-stream  Mach  number  of  6.2.  This 
comparison  is  presented  as  a  wall  temperature-time  history 
at  the  stagnation  point.  Since  the  test  model  required 
0.2  second  to  overcome  the  initial  temperature  lag  due  to 
the  starting  transients,  the  theoretical  analyses  per¬ 
formed  for  this  study  were  started  at  0.2  second,  using 
these  temperatures  as  the  starting  conditions.  As  seen 
from  figure  18,  the  theoretical  and  wind-tunnel  data  are 
in  good  agreement.  The  slight  deviation,  it  was  felt, 
was  due  to  either  inaccurate  temperature  measurements  in 
the  test  or  to  incorrect  skin  thickness  data  at  the 
thermocouple  locations. 

2.  Recommendations 

The  results  of  this  investigation  indicate  that  the 
heat  transfer  and  temperature  profiles  can  be  reasonably 
predicted  with  the  methods  used  in  this  study.  However, 
some  areas  could  use  further  study.  These  areas  are 
bluntness  effects  on  boundary-layer  edge  conditions  and 
heat  transfer  as  well  as  flow  conditions  in  the  region 
immediately  behind  the  stagnation  point.  Also,  the  ef¬ 
fects  of  multiple  shock  interaction  on  heat  transfer  in 
these  areas  are  a  possible  problem  for  study. 
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APPENDIX  A 


REYNOLDS'  ANALOGY 

(3) 

In  1874,  Reynolds  postulated  that  the  laws  govern¬ 
ing  momentum  and  heat  transfer  were  the  same.  In  1883  he 
discovered  the  laws  of  friction  resistance  for  turbulent 
flow.  He  was  able  to  express  the  analogy  between  heat 

transfer  and  momentum  transfer  in  mathematical  form. 

The  following  is  the  Reynolds  equation  of  analogy  to 
find  the  rate  of  heat  transfer  in  forced  convection  tur¬ 
bulent  flow  over  a  flat  plate.  At  a  point  along  the  plate 
X  distance  from  the  leading  edge,  the  equation  can  be 
written 

dT 

A  =  -T  C  t —  ( A— 1 ) 

Mcx  cx  p  du  1  ; 

By  rearrangement 

rua,  f  ^oo 

I  ^cx  du  =  -  L  TcxcpdT  (A'2) 

o  S 

After  integration,  it  becomes 

■*CXU-  -  tcxcp(TS  -  T->  <A-3> 

The  turbulent  friction  coefficient  is  defined  as  the 

ratio  of  the  combined  friction  force  x  to  the  total  ki- 

c 

2 

netic  energy  (1/2) pu  .  Therefore,  for  any  station  x, 
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and 


(A-4) 


4cx  -  hcx  (TS  -  T»>  <A-5» 

Substituting  into  equation  (A-3)  and  reducing 

hcx(TS  -  T«)u»  =  cfx  I  »u»cp<Ts  -  T»)  (A-6) 

or  by  definition  the  Stanton  number  is 


hcx  _  1 
Pcpu»  1  fx 


(A-7) 


It  must  be  understood  that  equation  (A-7)  is  valid  only 
for  Pr  =  1.  Thus 


(A-8) 


or 


CfRe 


(A-9) 


The  reference  temperature  is  used  for  evaluating  the 
physical  gas  properties.  Since  Reynolds  did  not  consider 
the  laminar  sublayer  in  his  analogy,  the  equation  is  valid 
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only  under  the  assumption  that  heat  resistance  through 

the  thin  layer  is  negligible. 

For  values  of  Prandtl  number  of  0.6  <P  <1.0, 

r 

Colburn, ^  to  obtain  more  accurate  data,  gave  the  follow¬ 
ing  equations  for  heat  transfer  in  turbulent  flow  over  a 

5  7 

flat  plate  with  5  x  10  <  Rq  <  10  in  terms  of  his  j-factor 


N  p2/^3  = 
ST  r 


wf 

T~ 


=  3 


(A-10) 


which  is  based  upon  the  friction  coefficient 

C*  =  0.0584  fT1/5 
r  e 


(A-ll) 


Combining  equations  (A-10)  and  (A-ll)  and  substituting 
into  equation  (A-8)  solving  for  the  Nusselt  number 


0.0296 


pl/3R4/5 
r  e 


(A-12) 


Therefore,  from  equation  (A-12),  the  local  convective  heat 
transfer  coefficient  for  turbulent  flow  over  a  flat  plate 
is 


hcx  =  °-0296  £  Re'8pr/3  (A-13) 

This  solution  is  applicable  to  conical  flow  when 

modified  by  the  Mangier  transformation^  and  for  compres 

12) 

sibility  by  Eckert's  reference  enthalpy  method.  So 
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for  conical  flow 


is  the  ratio  of  the  friction  factor  in  compressible  flow 
to  thd  friction  for  a  constant  property  value  gas. 


APPENDIX  B 


DETERMINATION  OF  BOUNDARY -LAYER 
EDGE  PROPERTIES 

To  calculate  the  convective  heat  transfer  rates, 
several  boundary-layer  edge  conditions  must  be  known, 
namely,  temperature  and  velocity  to  calculate  the  refer¬ 
ence  enthalpy,  h*,  and  pressure  to  calculate  the  refer¬ 
ence  density.  In  this  study,  these  parameters  were  ob¬ 
tained  by  curve-fitting  the  conical  flow  results  by 
Bertram and  references  17  and  13.  These  results  have 
been  correlated  as  a  function  of  the  hypersonic  similar¬ 
ity  parameter,  M  sin0  =  K  ,  from  which  the  following 
relations  were  obtained: 


The  local 


1 


1.4 


0.5 


2.5  +  8  K  ' 

_ c 

1  +  16  K 

C 

1  +  0.0966  K  +  0.2267 (K  )2 
c  c 

Mach  number  can  be  computed  from 


1  +  2.8  K2 
c 


V. 


/y  R  T 
e  e  e 


•^ooVoo 


(B-l) 


(B-2) 


(B-3 ) 


(B— 4 ) 
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From  previous  assumptions  and  reducing  equation  (A-4), 
the  following  ratio  is  computed: 


The  reference  enthalpy  is  obtained  from  the  curve 
fit  relation 


(B— 6 ) 


Assuming  y  =  1.4,  r  =  0.85  and  c  =  c  =  c 

Pg  Poo  Pyj 

jp  *  0.5  +  0.  +  0.0388  M*  (B-7) 

For  evaluating  the  fluid  properties,  the  reference 

temperature  is  used  which  is  expressed  as  a  function  of 

the  local  Mach  number  and  of  the  ratio  T^/Tg.  However, 
(18) 

Banner  et  al.  shows  that  the  method  of  Eckert  over¬ 

estimates  the  measured  levels  of  the  turbulent  heat  trans¬ 
fer.  They  suggest  that  this  conservatism  can  be  reduced 
by  neglecting  the  effect  of  heating  rate  in  the  calcula¬ 
tion  of  the  heat  transfer  coefficient.  This  is  accom¬ 
plished  by  substituting  the  boundary-layer  recovery  tem¬ 
perature  for  the  skin  temperature  in  the  equation  used  to 
calculate  the  reference  temperature.  The  equation  for 
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this  method  is  given  by 


TIw  =  l(Te  +  Tr  +  °*22  ^l^e)  (B'8) 

which  is  called  Eckert’s  Adiabatic  Wall  Reference  Temper¬ 
ature  Method. 


43 


r 


APPENDIX  C 

THEORY  AND  PROOF  OF  THE  FINITE  DIFFERENCE 
SOLUTION  FOR  FOURIER  CONDUCTION  EQUATION 

To  explain  the  theory  of  the  finite  difference  method, 
a  wall  is  divided  into  a  number  of  lamina  with  known  equal 
thickness,  Ax,  specific  heat,  and  thermal  conductivity. 

The  temperature  profile  within  the  wall  at  any  time  t  is 
represented  by  the  line  Tn_2>t  Tn.1<t  Tn<t  Tn+1(t  -rn+2ft- 
The  subscripts  n,  n  +  1,  etc.,  indicate  the  location  of 
the  points,  while  the  subscript  t  indicates  time.  For 
example,  the  temperature  of  the  wall  at  the  surface  n  in 
time  t  is  T  . .  See  figure  19.  At  a  time  interval  At. 
later,  the  temperature  at  the  same  surface  is  T  .  .  Ac- 

II  f  tT  1 

cording  to  this  method,  the  temperature  can  be  determined 
at  any  future  time  for  any  surface  location. 


Figure  19.  Schmidt  Plot  in  an  Infinitely 

Thick  Wall 
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To  prove  the  finite  difference  theory,  express  the 
finite  increments  of  temperature  in  the  Fourier  conduc¬ 
tion  equation  in  terms  of  temperatures  as  follows: 


A,  T  =  T  .  ..  -  T 
t  n,t+l  n,t 

(C-l) 

(AT)  =  T  ,  -  T  , 

x  ,  n,t  n-l,t 

(C-2) 

n-1 

(AT)  =  T  ,  .  ,  -  T 

x  'n+1  n+1 , t  n, t 

(C-3) 

A*T  =  (AT)  -  (AT) 
x  x  n+1  x  n-1 


*  <Vl,t  -  Tn,t>  -  (Tn,t  -  Vl,t> 


Tn+l,t  ^Tn,t  +  Tn-l,t 


(C-4) 


The  substitution  of  equations  (C-l)  and  (C-4)  into  the 
Fourier  conduction  equation  gives 


—  (T 

At v  n,t+l 


-  T  )  = 
n ,  t 


(Ax) 


2(Tn+l,t  -  2Tn,t  +  Vl,t>  (C-5) 


which  may  be  transformed  into 


2  (T  ...  -  T  ,  )  =  2a 
n,t+l  n,t 


At 


(Ax) 


T(Tn+l ,  t 


2T  +  T  ,  )  (C  —  6 ) 

n  / 1  n  x  / 1 


The  thickness  Ax  of  each  layer  of  the  slab  and  the  time 
increment  At  are  so  chosen 
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Then  equation  (C-6)  becomes 


(C-7) 


or 


Tn,t+1  -  +  Tn-l,t>  <c-8 

This  means  that  if  the  temperature  is  known  at  two  sur¬ 
faces  (n  -  1)  and  (n  +  1) ,  at  a  time  t,  the  temperature 
at  surface  n  midway  between  (n  -  1)  and  (n  +  1)  and  at  a 
time  (t  +  At)  is  equal  to  the  mean  value  of  Tn+1  fc  and 

T  ...  Since  the  intersection  of  surface  n  with  the 
n-i,  t 

line  joining  Tr+1  and  Tn-1  gives  the  mean  value  of 

these  temperatures,  the  theory  of  this  method  is  there¬ 
fore  proved. 
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APPENDIX  D 

PROGRAM  LISTING 
SLEAT 
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PROGRAM  SLEAT1TNPUT, OUTPUT) 

C  JAVELIN- J A VFL IN  N1TH  AEROOtNAHIC  HEATING 
COMMON  /A/  XL<?) «MAREA(3) ,AREA<3) 

common  /b/  5mich,RE 

CMflQN.  /c/  THETC.THETH.V.TEMP.AHACH.P 
COMMON  /O/  TAN,$TAN,OELT,RHO,VIS 

COMMON  /E/  TN(4,5)  .HUM  ,STM(54)  ,Q00TC<4)  ,SQDOTC ,RNC5)  ,NAL  ,TOL ,NAL 
l$,NOFX,STG,VIS€,VISG 
COMMON  /H/  AMACHE,PE,VE,TEMPE,RHOE 
PRINT  1*5 

PlHa.lW9—  - 

THETC=9.0333*P 1/180. 


NOFX-4 

NOFX«*»NOF**1 

NAL=3 

HALS* 54  .  _ 

TOL-.0104 

DATA  HARE  A/30 • 055/ 

OATA  XL/12.756/ 

DATA  AREA/, 44177/ 

DATA  RN/2.E6,2.E6,5.E5,5.E5,1.5E5/ 

. .  QAtA.  HL/2.»_jL»  ,.»5B11«»2IL/ 

TEMP-534. 

TAH=TEHP 
STAH*TEMP 
00  t  T*1,N0FX 
00  1  J-ltNAL 

1  T  til  J  «1L*  T.EMP- . . . . 

00  2  I=1,NALS 

2  STH(I)*TEMP 

VISa7. 285E-0  7* (TENP**1. 5) / (TEMP*198. 72) 

N5TG*1 

GRAVsJlZ^A. . . . 

GF*0. 

X=0.0 

TIME=0.0 

V*0.0 

JJ*0 

P~  25  *.75527.1*13  . . . . 

R*53. 3 

RHO=P/(R*GRAV*TEMP) 

VS=49.06*SQRT(TEMP> 

TRST=0. 

ORAG-O . 

c.br?o,q.  ...  .  _  ..  .  . . 

co=o. 

COH=0, 

CF*0. 

COT0T=0. 

VMAS=16.687l 

PRINT  16 _ _ _ _ _ 

PRINT  17,  TTNE,QRAG,TRST,VMAS,V,X,ANACH,GF,CBR,CD,CON,CF,COTOT 
PRINT  16,  f (TH(I,J),J*1,N0FX> ,1-1, NAL) 

PRINT  19,  <STMI>  ,I*1,NALS) 

PRINT  20,  COOCTC(I) ,I=1,N0FX) ,SQOOTC 
OFLTs.Ol 

L  ....  -TLMJEfTlHEtJSELL _ _ _ 

AHACMsV/VS 

JJ=JJ*1 

IF  <TIHE-l,95)  8,9,9 
21  IF  CTTME-1.85)  8,4,4 
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NSTG=  1 

IF  CTIHE-1.86)  5,6,6  _ 

NTL*3.1B2 
GO  TO  7 

6  NTL*3,9512#0ELT 

.  _T*  5 _ _ _  _ 

CALL  TPSTJ  CTRST,T> 

GO  TO  10 
T*TIME 

CALL  TPSTJ  CTRST,T) 

WTL--4.14469*0ELT 

__  GO  TO  JJ3 _ 

0  TRST*0. 

VHAS*9. 01973 
NTL=0 • 

CALL  COBRAKE  (C0R, AHACH) 

CBR*CBP*. 175/AREA 

COB*l.  125/AHACH**2 _ 

IF  C AHACH. LE .1.01  COB*  t, 125 

10  VHAS*VMAS~NTL 

IF  IAHACH.LT. .1)  GO  TO  U 
CALL  COSKIN  (CF , NSTG, V ,RHO,VIS) 

11  CALL  COHACH  ( AHACH, CO ,CDN, V, NSTG) 

_  COTOT*  CC0»C0N»C08»CF»CHR) _ _ 

ORAG*.5*RHO*AREACNSTG)*COTOT*V**2 

NT*VMAS*GRAV 

GF*CTRST-ORAG)W 

V*VMTRST-ORAG)*OELT/VMAS 

X*X»V*OELT 

_ IF  (AHACH.LT, l.Q)  GO  TO  1? _ 

CALL  HEAT 

12  IFCJJ.LT. 100)  GO  TO  13 

JJ*0 

PRINT  17,  TIHE, DRAG,  TRST , VMAS *  VfXtAHACHfGF,CBR, CO, CON, CF, COTOT, COB 
lySTGyPEyAHACHE, VIS, VISE ,VJSG 

PRINT  18.  MTMI.J)  ,J«l,NOFX)  .  1*1. NAL) _  _ 

PRINT  19,  CSTNCI) , 1*1 , NALSl 

PRINT  20,  (QOCTCCI)  ,I*l,NOFX) , SQOOTC 

13  IF  (V.LE.O.)  GC  TO  14 
IF  CTIHE-5.)  3,3,14 

14  CONTINUE 

_ C _ _ 

15  FORMAT  C1H1I 

If  FORMAT  COSH  TIHE  DRAG  TRST  VHAS  V  X 

1  HACH  GF  CBR/42H  CO  CON 

2  CF  COTOT  ) 

17  FORMAT  C//F10.2,8E14.4/10X,9E14,4/10X,9E14.4) 

_ 18_  FORMAT  C19H _ TN  ,4F8 . 0/C  19X.4FB.  0)  ) 

19  FORMAT  C23H  TNCSTAG)  ,1SF6.0/(23X,18F6,0~>  > 

20  FORMAT  C34H  QOOTC  CNOFX) , SQOOTC  ,5E14.4) 

ENO 


l 
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SUBROUTINE  HFAT 
C  CONE  HEATING 

COMMON  tZf  thet^,theth,v,teh?,ah,? 

COHNON  /O/  TA¥»  ST AH ♦ CELT  »RHO*  VIS 

COMMON  /T/  TW(4,G),WLT4),STWr54T,(rD0TC(4),50D0TC,RN(5>,NAL,T0l,NAL 

1 S  t  NOF  t ,  GTG>  VISE  »  V I  S  L, - 

OIMENSION  VRHO(4,5),  VTC(4,5),  VCP(4,5>,  SACP(54),  SATC(54>,  SAPHO 
t  (541 
AMAOHs AM 

CALL  BOLAT  (TEMPE* °E» VE •  ANACHE  tRHOE) 

PTS=P*  (  (1.2*AMACH**?)**3.51»U6.a/  ( 7 . 0 • AHACH**2-1 . 0) >  **2.5) 

-  -tjtr  t~l  FW-i ,  vny  t - 

KM 

TR=0. 5*  (TEMPF*TPfC»n.4  4*RE*AHACHE**?»TFHPE) 

CALL  TCAT  (TCA.TR) 

CALL  VISCOT  (VTSG.TR) 

CALL  VISCOT  (VISF,TEMPE) 

- PWT=8J.98T - - - 

GRAVM2.174 

GsGRAV 

AJ=779. 

CPA=.27P5 

rhotmits. 

- gHOA^lftO. - 

RX=.341f 

SX=.0ir45*RX*PMI 
PP=CPA*VISG/TCA 
MNRE*RMOf *VE*NL  (LEN)  m?E*G 
RF  *PR* • . 33 

- tttos»i»rr  trrtfr;  »t*) - 

HE-.3  33SMCA/NL  (LEN)*PR**.yj*HNRE*T.8*  (VISG/VISE)  (RHOS/RHOE)  * 

1  ••  8 

OTR=VE*M/(?.*G*AJ*CPA) 

TAM=TEm»E»RF*nTP 
TREC*TENPE»RE* (STG-TEMPE) 

- QOOTC  tLFNI  =Wf  f  TAWTWTr.Lfim - 

1  THKsTHTK, LENT 

IF  aEN.GF.NOFTO  DO  TO  ? 

CALL  ALUMTC  <TC  * TWK) 

call  aluhcp  (CP) 

VRHO(K,LEN) *PH0A 

- en  T0  ^ - 

?  CALL  TUNGTC  (TC,TH<) 

CALL  TUNG CP  (C«»,TWK) 

VPHO(K,LPN>=PHOT 

3  VTC (K  t  LEN) -TC 
VCP(K,LEN)=CP 

- Ty  BO  TO  4 - 

IF  (K.FO.NAL)  GO  TO  6 

P1=VTC(K-1  ,LFN)  *  (TH(K-l,LEN)-TH(K,LEN)  ) / (VRHO (K-l,LFN)  *VCP (K-1,LFN 
t ) • TOL •  •?) 

P7-VT3  ( K  ,  LFN)  *  fTW(<,LEN)-TH(K*IiLEN)  )  '  tVRHO  (K  ,  LFN)  *VCP(KfLEK)  *TOL* 
1*?) 

- nrfr,Lrir>  gTwnertFNi  aoelt*  (pi-p?> - 

GO  TO  3 

4  P3  =  0ELT/(VCP(KtLEN)  *VRHC(K,LENJ  *TOL)*2. 

P4=HE* (TAH-TW(K,LEN) ) 

?3=VTC  CK,LFN)  MTWt.(,LFN)-TNTXM,LEN)  )  FTCL 
TW  (K»  LEN)  =TWOC,LEN)  APT*  (P4-P5) 

- , - mcrt - 

GO  TO  1 

6  TH(RtLEN)  =TWTK,LFN)  +OELT*  (VTC ( K- I t LEN)  /(VRH0(K-1  f  LEN)  *VCP(K-1  ,LEN> 
l*TOL*»?)MTH(K-l,LEN)-TH(K,LEN))  ) 


r 


7  CONTINUE 
KS  =  1 

LEN=NOFXM 

SDTR=V*»?/(?.*G*AJ*CPAI 

STG  =T  SDTR  - - 

- ST^5TEHPf,5S*lS!Hlll-!  kBrJM-.iy*  rSTG-lkPHI - 

CALL  TCAT  (STC A » STR) 

CALL  VTSCOT  (SVISG.STR) 

SPR=CPA*SVTSG/STCA 

SRF=SPR**.T! 

RH0WS=PTS/(1716.*STR) 

- HS  =  Tk  HP  »  OTTTTT 2 YTE^TE.1* TTWPrTT\TrT?7~ry7rT7yTJ7 - 

HW=STW(l)*(0.23+lZ.5E-OE*STH(in 

GTAR=(  (1.0  +  5. 0/(AMACH**?))  *(1.0-1. 0/(1.  4*AMACH**2>  >  >**0.25 
HSE=V*»?/ (2.*G*AJ) 

S0DnTC=0.94*(RH0HS*SVISG*V*G)**.5*HSE*GTA8/(RX**.5> 
SHF=SQnOTC/(STG-5TW(l)  ) 

- ST3R-_-tEWPI.w,t.u^ - 

STRFC  =  TFMP+SRF»  (STG-TEMP) 

8  STHK=STW ( KS) 

CALL  TUNGTC  TTvT »TTHX7 
CALL  TUNGCP  (C",STWK) 

SARHO  (KS)  =RHfT, 

- yfl-CP  TXS'  T  J:UP - 

SATC(KS) sTC 
IF  (KS.E0.1)  GO  TO  9 
IF  (KS • FO.NALS)  GO  TO  11 

STWUS)=STW(KS)*DELT*TSATCTK5-i)  / IS  ARHO  (KS-1)  *SACP  t KS-1  >  *TOL**2>  •  ( 
lSTW(KS-l)  -STW(KS))-SATCfKS)y(S*RHO(KS)*SACP(KS>*TOL**2>*(STW(KS>-S 

- 7TH(KSUTn - 

GO  TO  13 

9  STH (KS)  =STW  (KS)  VOELT/rSSITPTKST  ^S  ARHO  (KSI »TCL)  *  (SHE*  (STG-STW  (KS)  )-S 
1ATC  (KS)  /TnL*  (STWTK^-STWTirSM)  )  I  *2. 

1)  KS=KS*1  “ 

GO  TO  8 

n - 5TH(KS1-,STHTKS) THE  STTC  TXS-1 )  /  ( SA'RHO'(KS-l)*  STCP  (KS- IT  *TDL*  *?)'*! 

ISTW(KS-  1) -STWrK^ll  - 

RETURN 

ENO 
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APPENDIX  F 

HEAT  TRANSFER  TO  WEDGES 

Heat  transfer  In  special  areas  such  as  in  the  regions  of  shock-wave  boundary 
layer  interactions  and  in  notches  or  steps  is  extremely  difficult  to  analyze 
theoretically.  Therefore,  some  simplified  relationships  for  approximating  the 
heat  transfer  in  these  areas  would  seem  to  be  in  order.  The  relationships  pre¬ 
sented  in  this  appendix  for  the  heat  transfer  coefficient  are  intended  to  only 
serve  as  guides,  since  for  the  most  part,  verifying  experimental  data  are  lacking. 
The  heat  transfer  to  a  surface  in  the  vicinity  of  a  notch  or  step  is  influenced 
by  the  change  in  flow  field  caused  by  the  notch  or  step  and  heat  transfer  coeffi¬ 
cients  tend  to  increase  in  these  regions.  These  relationships  are  somewhat  con¬ 
servative;  however,  use  of  them  in  localized  areas  will  not  influence  overall 
design  weights  significantly. 

The  heat  transfer  coefficient  for  the  surface  of  an  upstream  or  forward  facing 
protuberance  can  be  calculated  from 

h  =  0.0296  -  P  •  33  R  -8  (F-l) 

e  x  r  e 

where  the  Prandtl  number,  Pf,  and  the  Reynolds  number,  R^ ,  are  evaluated  at  local 
boundary  layer  conditions  for  the  same  location  with  no  step.  The  thermal  con¬ 
ductivity,  k,  of  air  is  evaluated  at  the  reference  temperature. 

The  wedge  is  also  divided  into  a  number  of  segments  similar  to  the  cone,  but 
the  wedge  is  assumed  to  be  a  flat  plate  inclined  to  the  stream  flow.  Since 
there  is  an  interaction  of  the  nose  cone  shock  wave  with  the  wedge  and  the  wedge 
is  completely  within  the  shock  wave,  it  was  assumed  that  the  ambient  conditions 
for  the  wedge  are  the  cone  boundary  layer  properties.  In  other  words,  the  wedge 
is  assumed  to  be  exposed  to  the  heated  environment  corresponding  to  the  cone 
boundary  layer  edge  conditions. 
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